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The soybean cyst nematode 
(SCN), Heterodera glycines, is 
the most damaging pathogen 
of soybeans in the United 
States. To assess the severity 
of nematode infestations in 
the field, SCN egg population 
densities are determined. 
Cysts (dead females) of the 
nematode must be extracted 
from soil samples and then 
ground to extract the eggs 
within. Sucrose centrifugation 
commonly is used to separate 
debris from suspensions of 
extracted nema- tode eggs. 
We present a method using 
OptiPrep as a density gradient 
medium with improved 
separation and recovery of 
extracted eggs compared to 
the sucrose centrifugation 
technique. Also, computerized 
methods were developed to 
automate the identification 
and counting of nematode 
eggs from the processed 
samples. In one approach, a 
high-resolution scanner was 
used to take static images of 
extracted eggs and debris on 
filter papers, and a deep 
learning network was trained 
to identify and count the eggs 
among the debris. In the 
second approach, a lensless 
imaging setup was developed 
using off-the-shelf 
components, and the 
processed egg samples were 
passed through a microfluidic 
flow chip made from double- 
sided adhesive tape. 
Holographic videos were 
recorded of the passing eggs 
and debris, and the videos 
were reconstructed and 
processed by custom software 
program to obtain egg counts. 
The performance of the 


software programs for egg counting was charac- terized with SCN-infested soil 
collected from two farms, and the results using these methods were compared 
with those obtained through manual counting. 


Introduction 


Many nematodes (microscopic roundworms) are soil-dwelling plant 
parasites that infect the roots of plants and cause billions of dollars of 
crop loss worldwide on an annual basis. Hetero- dera glycines, the 
soybean cyst nematode (SCN), is the most damaging soybean pathogen 
in the United States and Canada, causing hundreds of millions of dollars 
in crop loss annually [1]. The SCN infects the roots of soybeans and 
siphons nutrients from the plants, leading to stunted 
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growth and reduced crop yields. The amount of damage and yield loss 
caused by the nematode is related to several factors, including the 
number of eggs (egg population density) in the soil. Knowing the SCN 
population density in the soil can be useful in guiding the use of 
management strategies and assessing the success of management efforts 
in SCN-infested fields. 

To determine the population density of SCN in a field, one or more 
multiple-core soil sam- ples are collected. On a sample-by-sample basis, 
the cysts (egg-filled dead SCN females) are extracted from samples, and 
the eggs are extracted from the cysts, then counted using a micro- scope. 
Two methods commonly used to extract nematode cysts from the soil are 
(i) wet sieving and decanting, and (ii) elutriation. Wet sieving and 
decanting [2] involve suspending soil in water, agitating the suspension, 
allowing the heavier soil particles to settle to the bottom of the container, 
and pouring the suspension through two sieves. The top sieve (usually 
with 850- um-diameter pores) will capture root fragments and other 
debris, which will be discarded, and the bottom sieve (usually with 250- 
um-diameter pores) will capture nematode cysts and cyst- 
sized debris. Smaller debris (<250 um in diameter) will pass through both 
sieves and be dis- 
carded. For the elutriation method, soil is suspended in a column or cone 
of upward-flowing water. Heavy soil particles remain near the bottom of 
the flowing water suspension, while cysts and less dense objects float and 
pour out of the top of the column or cone to be captured on sieves [3] as 
with the wet-sieving and decanting method. 

Currently, two methods commonly used for extracting eggs from 
nematode cysts are (i) by grinding the cysts in a plastic or glass tube with 
a stainless-steel pestle [4] or Teflon tissue homogenizer and (ii) by 
grinding the cysts on a 250-um-pore sieve with a rubber stopper [5]. Both 
methods result in the capture of eggs and similarly sized debris on 25-um- 
pore sieves. 

Accurately and efficiently counting the eggs to determine the population 
density can be diffi- cult and inefficient when considerable amounts of 
debris are recovered with the eggs. Separat- ing the eggs from the debris 
often is necessary to make counting possible. Also, it is desirable to 
separate and discard debris from eggs in suspension when the eggs are 
to be used in labora- tory experiments, such as when assessing 
development and studying hatching of nematode juveniles from the eggs 
[6-8]. 

Sucrose centrifugation [9] is a method commonly used to separate 
plant-parasitic nematode juveniles and eggs in suspension from debris. 
Although sucrose centrifugation is inexpensive and easy, it may not be the 
most efficient method and it may have adverse effects on the nema- todes. 
Exposing nematode eggs and juveniles to sucrose solutions, with a high 
osmotic potential, may harm the eggs and juveniles if exposure to the 
sucrose is prolonged. Also, if the nematodes are not thoroughly rinsed 
with water following centrifugation in sucrose solution, sucrose resi- dues 
may remain and promote bacterial and fungal growth on the eggs and 
juveniles. Deng 
et al. [10] suggested an alternative density gradient method to sucrose 
centrifugation using an iodixanol solution, OptiPrep” (also known as 
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Visipague™ in medical uses). Testing the effects of a single concentration of 
OptiPrep", they found that the efficiency of extraction of the reniform 
nematode (Rotylenchulus reniformis) and post-extraction mobility of the 
recovered nematodes were both 100% greater than when the sucrose 
centrifugation method was used. 

Following the extraction and cleaning of egg samples, the eggs are 
counted to determine the population density of nematodes. Counting is 
completed manually through microscopic obser- vation, requiring 
trained personnel. The process is time intensive, laborious, and prone to 
human error. It would be ideal if this portion of the process could be 
automated to address these shortcomings. Furthermore, automated egg 
counting could reduce labor costs and processing fees compared to 
manual counting methods where the expense increases with sample 
quantity, 

The objectives of this work were: 1) to develop a method of purifying 
eggs from suspension with debris and 2) to develop new methods to 
automate counting of eggs once extracted from cysts and purified. In the 
first counting approach, a high-resolution scanner takes images of 
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the processed sample (i.e. stained eggs with debris) dispersed on filter 
paper that was then run through deep learning algorithms to 
automatically identify and count the eggs. In a second counting approach, 
a benchtop, lensless imaging setup takes real-time, holographic videos of 
the processed sample passing through a flow chip and is analyzed with a 
custom software pro- gram to determine the egg count. We performed 
detailed characterization of the new methods, while attempting to 
minimize the number of manual steps. 


Materials and methods 
Sample preparation 


Soil samples were collected from two fields in Muscatine and Story 
County in Iowa. The fields were located on research farms owned and 
operated by Iowa State University. No specific per- missions were 
required to obtain the soil because the purpose of the farms and their 
fields is to support crop research, and collection of soil samples is a 
routine part of conducting such research. The field studies did not involve 
endangered or protected species. The sample prepa- ration involved two 
steps: egg extraction and egg staining. 

To extract the eggs, each soil sample was poured in a bucket filled 
with water (approxi- mately 2 liters), mixed thoroughly, and allowed to 
settle. The soil suspension then was poured through a 20-cm-diameter 
sieve with 850-um-diameter pores above a 20-cm-diameter sieve with 
250-um-diameter pores [2]. The debris and egg-filled cysts (dead 
nematode females) were collected on the 250-um-pore sieve and then 
transferred to a 3.7-cm-diameter, 250-um- pore sieve and were crushed 
using a motorized rubber stopper to release the eggs [5]. These eggs 
(along with similar sized debris) were collected on a 15-cm-diameter 
sieve with 37-um- diameter pores and then transferred into a 
microwavable container [11]. 

For staining the eggs, a stain solution was prepared by adding 3.5 g 
acid fuchsin (F8129, Sigma Aldrich) and 250 mL glacial acetic acid 
(ARK2183, Sigma Aldrich) to 750 mL of distilled water and stirred well [4]. 
One drop of the stain and four to five drops of 1 M HCl solution were 
added to each plastic beaker containing the eggs. The beakers and their 
contents were heated in a microwave for 15 seconds to stain the eggs 
[12]. 


Sample cleaning method using OptiPrep™ -based 
centrifugation 


Fig 1A is a schematic representation of the protocol for sample cleaning 
using OptiPrep" as the density gradient medium. The protocol consists of 
two steps: centrifugation and separa- tion. During the centrifugation 
step, 3 mL of the density gradient medium (OptiPrep”) at a specific 
volume percentage was put in a 15 mL centrifuge tube (Corning"). 
Thereafter, 5 mL of the stained egg sample was pipetted on top of the 
OptiPrep" solution, forming an emulsion interface layer between the 
solutions of different densities as shown in Fig 1A(i.). The com- bined 
sample was centrifuged at 840 G for 2 minutes. Particles with higher 
density than the gradient solution passed through the emulsion layer and 
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were deposited in a pellet at the bot- tom. Particles with comparable 
density to the gradient solution were concentrated in the inter- face 
layer, while particles having lower density floated in the top layer of the 
solution. The contents of the centrifuge tube were categorized into three 
layers: top (4 mL), interface (2 mL), and bottom (2 mL). In the separation 
step, each of the three layers were pipetted into three dif- ferent 15 mL 
test tubes and diluted to 11 mL with water, as shown in Fig 1A(ii.). 


Scanner-based egg counting method 


Hardware and software components. Fig 1B illustrates the scanner- 
based egg counting method. Here the three distinct layers of processed 
sample (top, interface, and bottom) obtained 
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Fig 1. Overview of the new methods of sample cleaning and egg counting. a) Sample cleaning used a density gradient 
medium (OptiPrep") and consisted of two steps: centrifugation and separation. i) An illustration of the density gradient 
centrifugation is shown where the three distinct layers (top, interface, and bottom) were visible after the centrifugation step. ii) The 
three layers were separated and diluted with water to 11 mL. b) A scanner-based egg counting method is shown where three 
solutions were dispersed on separate filter papers, allowed to air dry, and then placed face-down on an A4-sized (21.6 x 27.9 cm) 
transparency sheet. The filter papers were scanned at 4800 dpi and analyzed by a deep learning network model to obtain a count of 
the eggs. c) The lensless setup consisted of a light-emitting diode (LED) to illuminate the processed sample flowing through a 


microfluidic flow chip and a CMOS image sensor to record videos of the eggs in real-time. A custom software program, written in 
Python, analyzed the videos and produced the egg count. 
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after centrifugation in OptiPrep” were poured on separate circular filter 
papers (Grade 41, What- man”, 90 mm diameter) and allowed to air dry 
at room temperature (approximately 22°C). The filter papers were placed 
on an A4-sized (21.6 x 27.9 cm) transparency sheet face-down and placed 
on a flatbed scanner (Epson Perfection v750 Pro) connected to a desktop 
computer. Using the scanner software (Epson Scan software, Version 
3.921), each piece of filter paper was selected with a bounding box and 
scanned as shown in Fig 1B. The following settings were chosen: reflec- 
tive scanning mode, 4800 dpi, and 24-bit color. All other imaging 
parameters were set to their default values. The scanned images were 
saved as JPEG files to the computer hard drive. 

Data collection and pre-processing. Each scanned image of the 

filter paper 
(17759 x 17759 pixels) was split into 4900 patch images (256 x 256 
pixels). The colorspace of each patch image was converted from RGB 
(Red, Green and Blue) to HSV (Hue, Saturation, and Value) to distinguish 
the different colors because the HSV colorspace separated the color 
information (chroma) from the image intensity (luma). A range of HSV 
values was selected to identify objects of a similar color to the stained 
eggs within the patch image. Then by applying thresholds on the object’s 
physical dimensions (i.e. width, height, shape, and area), all of the eggs 
were detected in the patch image. A label image with the eggs was 
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created and stored along with the patch image on the hard drive. The 
above process was repeated for all the patch images, and subsequently 
used for training and testing the deep learning network. The patch and 
label images were resized to 128 x 128 pixels to reduce the number of 
training parameters and model size. The pixels were normalized to ensure 
that all of the features were given equal importance. The total data set 
consisted of 60 filter paper images from 20 different soil 
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samples, yielding 294,000 patch images, which were randomly divided 
into training and test data sets (80:20 split) to be used to develop the 
deep learning model. 

Deep learning network architecture. To automatically learn the 
features related to the SCN eggs from patch and label images, we have 
employed a convolutional autoencoder network 
—a specific type of autoencoder network that used convolutional layers 
[13,14]. In general, a ‘convolutional network’ is composed of four layers: 
convolution, activation, pooling, and drop- out layers. The convolution 
layer had neurons with weights and biases, which were updated after 
every iteration by a backpropagation algorithm in the training process. 
The activation layer con- sisted of a non-linearity function and performed 
mathematical operations on the input. An example activation function is 
the ReLU (Rectified Linear Unit) and has been shown to greatly 
accelerate the training process [13]. The pooling layer performed non- 
linear down-sampling on the input image by extracting the maximum or 
average of all the non-overlapping sub-regions in the image. This layer 
reduced the number of training parameters and memory footprint of the 
network. The dropout layer selected a random set of neurons determined 
by a percentage proba- bility and set their inputs to zero making them 
unusable in the decision-making process of the network. Besides the 
convolutional network, the ‘autoencoder neural network’ presented an 
unsupervised learning platform to provide an approximate mapping of 
the inputs and outputs. It consisted of the ‘encoder’ and ‘decoder’ paths. 
The ‘encoder’ path compressed the input informa- tion by downsampling 
and learned important features, while the ‘decoder’ path reconstructed an 
approximate higher dimensional output utilizing upsampling operations. 

A schematic of our network architecture is shown in Fig 2, which is a 
modified version of the U-Net convolutional autoencoder model [15]. We 
added dropout layers to prevent overfit- ting, replaced unpadded 
convolutions with padded convolutions to avoid cropping operations, and 
replaced the soft-max activation with sigmoid activation. We reduced the 
feature maps to half of the original size, which decreased the number of 
parameters. 

As shown in Fig 2, the deep learning network took 3-channel RGB 
patch images as inputs and produced 1-channel grayscale label images as 
outputs. The network consisted of nine layers; each layer had two 
convolutional layers and a max pooling layer. In the ‘encoder’ path, each 
step performed repeated convolutions with a 3x3 kernel followed by 
ReLU activation function and a 2X2 max pooling operation. In addition, 
the fourth and fifth layers contained dropout lay- ers with probabilities of 
0.4 and 0.3, respectively. The filters were doubled as the downsampling 
operation was performed to move to the next layer. The number of filters 
started from 32 in the first layer and extended up to 512 in the fifth layer. 
Thereafter, the ‘decoder’ path was initialized wherein, at each step, an 
upsampling operation and 2x2 convolution was performed and fol- lowed 
by concatenation of the corresponding feature map from the ‘encoder’ 
path. Then, simi- lar to the step in ‘encoder’ path, repeated convolutions 
with a 3x3 kernel followed by ReLU activation was performed. A final 
1x1 convolution and sigmoid activation was performed to map the label 
image. Details of our deep learning model are listed in S1 Appendix. 
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Training and testing the deep learning network model. The 
patch images and their cor- responding label images were used to train 
the network with Adaptive Moment Estimation (ADAM), which is a type 
of stochastic optimization of cross entropy loss [16]. In total, the net- 
work had 7,760,097 trainable parameters and was processed using an 
NVIDIA™ Geforce GTX 1070 Ti graphics processing unit (GPU) with 8 
gigabyte memory and 2432 CUDA cores. The batch size was set to 32 
and the training was performed for 100 epochs. The model imple- 
mentation was written in Python using Keras [17] and Tensorflow [18] 
libraries. 

The test dataset comprising the patch images was run through the 
above trained model and the corresponding grayscale images were 
generated. A thresholding technique yielded binary images which were 
passed through a blob labelling algorithm to get the egg counts in the 
patch images [19]. 
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Fig 2. The convolutional autoencoder network model. Each box represents a three dimensional feature map (W, H, and D). 
The depth (D) of each box is mentioned on top and the x-y sizes (W, H) on the bottom. The blue box represents the feature map 
copied from the encoder step and is concatenated with the feature map (black box) generated by upsampling the previous layer. 
The input to the network is a 3-channel RGB image and the output is a 1-channel grayscale image. The arrows represent different 
operations. 
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Lensless imaging method 


Hardware and software components. Fig 1C shows the setup for 
the lensless egg count- ing. The basic principle of lensless imaging is 
discussed elsewhere [20], and has been used to image microscopic 
objects such as cells, bacteria, and even nematodes without the need for 
expensive microscopes. Our lensless imaging hardware consisted of the 
lighting, imaging, and processing modules. (i) The lighting module 
comprised a light emitting diode (LED, wave- length A = 616 um, Vishay 
Intertechnology) that was aligned with a 100 um pinhole (Edmund 
Optics) and sealed within a pinhole mount. This arrangement allowed the 
light from LED to be emitted through the pinhole producing spatial 
coherence [20]. The LED was connected in series with a 60 QO resistor 
and is powered by the 5 V GPIO (general purpose input/output) pin on 
the Raspberry Pi (RPi) 3 Model B board. Fig 3A shows the wiring 
diagram. (ii) The imaging module comprised a Pi camera (8 MP) with a 
CMOS sensor. The default housing of the cam- era was removed using a 
standard razor blade to expose the CMOS sensor. Fig 3A shows the 
camera connected to the CSI (camera serial interface) port on RPi. (iii) 
The processing module comprised a portable RPi microcomputer having 
a 1.2 GHz 64 bit central processing unit (BCM2837), 1 GB random-access 
memory, and a built-in Wi-Fi unit. The RPi ran on the Raspbian Jessie 
operating system (Debian “Jessie” based/Linux kernel) loaded onto a 64 
GB microSD card (SanDisk"). Camera support was enabled in the 
settings, and open-source video streaming software was installed to 
access the live video feed remotely [21]. The RPi was pow- ered by a 
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micro USB wall-mount power supply (5 V, 2.1 A). A 3-D rendering of the 
platform that housed the various modules is shown in Fig 1C and the 
platform served to reduce interfer- ence from ambient light. The platform 
was printed on a 3-D printer (da Vinci 1.0 Pro, XYZ- printing) using a 
Polylactic acid filament (1.75 mm diameter, XYZprinting). 
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Fig 3. Wiring diagram for lensless imaging and design of the microfluidic flow chip. a) The CMOS sensor was connected to 
RPi through the CSI port. The 5V GPIO pin of RPi was connected to the anode of the LED while the ground pin from RPi was 
connected to the 60 Q resistor. The cathode of the LED was connected to the other end of 60 O resistor. b) (i) An illustration of the 
microfluidic flow chip is shown. The microfluidic flow design was cut on a piece of double-sided tape using a cutting machine, and 
the tape was bonded to a microscope slide. A transparency film with input and output accesses was cut and adhered to the top side 
of the tape. Two circular tape pieces with holes cut in the center were attached to the fluid ports and placed around the access 


holes in the transparency sheet. (ii) A diagram of the assembled microfluidic flow chip is shown with all of the layers bonded 
together. Scale bar = 20 mm. 
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After sample cleaning as described earlier, 1 mL of the processed 
sample was loaded in a 3 mL syringe (BD Biosciences) and connected to 
a syringe pump (KDS-100, KD Scientific). The standard needle of the 
syringe was replaced by a dispensing needle (16 gauge5.08 cm, Howard 
Electronic Instruments) to avoid clogging. The syringe was positioned 
vertically and taped with a vibratory motor (3 V/60 mA, 6500 RPM, 
Jameco Reliapro) to prevent settling of debris in the tubing. The 3.3 V 
GPIO pin was used to power the vibratory motor. A plastic tubing (inner 
diameter: 16 gauge) dispensed the processed sample from the syringe to 
a microfluidic flow chip. The RPi was turned on and video recording was 
enabled as the eggs and debris passed through the flow chip. Another 
piece of plastic tubing was connected to the output port to direct the 
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processed sample into a waste reservoir. After all of the processed 
sample passed through the flow chip, the video recording was 
terminated. 
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Design of the microfluidic flow chip. The microfluidic flow chip 
was intended to be made using low-cost and simple-to-use materials and 
tools, thereby eliminating the need for micromachining or 
microfabrication techniques such as photomask design, lithography, spin 
coating, developing, curing, and etching. The flow chip comprised three 
layers: the base, chan- nel, and cover layers. A microscope slide (25 x 75 
x 1.0 mm, Fisherbrand") served as the base layer. The channel layer was 
designed and cut from double-sided tape (25 x 75 x 0.2 mm, 3M"). The 
cover layer was constructed from a transparency sheet (Apollo"). Initially, 
commercial soft- ware (Studio", Silhouette America) was used to create 
the design of the microfluidic flow chan- nel (1.4 x 50 x 0.2 mm) witha 
central imaging chamber (2.1 x 2.4 x 0.2 mm). The different layers are 
shown in Fig 3B(i). Double-sided tape was attached to the cutting mat 
and loaded into the cutting machine (Cameo”", Silhouette America) to cut 
the flow channel onto the tape. 


The cut tape was then removed, aligned, and bonded to the base 
microscope slide. The design 


of the cover layer had holes for an inlet and an outlet which were aligned 
with the ends of the flow channel. The same cutting process described 
above was used to create the cover layer in a transparency sheet, which 
was placed and bonded to the channel layer. PDMS ports were used as 
input and output ports and sealed to the transparency film using circular 
rings of double- sided tape as shown in Fig 3B(i). The assembled flow 
chip is shown in Fig 3B(ii). 

Data collection and pre-processing. The videos recorded from the 
RPi were recon- structed using a custom software program written in 
Python. As shown in Fig 4, the program loaded one video at a time and 
read all frames sequentially. The objects of interest (eggs) in a frame 
were identified by subtracting the current frame from the previous frame 
and filtering the noise using a median filter. This step eliminated any 
static content in the current frame including imperfections in the flow 
channel, settled debris, and image noise. Only the moving objects 
remained after this step, including eggs and debris. The resultant image 
was in a binary format with object pixels colored in white against a black 
background. From the binary image, different physical dimensions (i.e. 
width, height, shape, and area) of the identified objects were determined 
and compared with those of a typical SCN egg. This helped identify if the 
selected object was a SCN egg. The process continued until all of the 
detected objects in the image were validated. The program then read the 
next frame and repeated the process until the final frame was reached. 
Thereafter, the total SCN egg count was reported for the processed 
sample. 

Principle of reconstructing holographic images. The spatially 
incoherent light from the LED passed through an aperture resulting in a 
partially coherent light called the reference wave U(x, y, Z). The 
reference wave illuminated the objects located at distance Zz, in the 
channel of the flow chip. The incident light was scattered to generate the 
object wave U(x, y, Z). The interference between the reference wave and 
object wave produced holographic patterns. The intensity /(x, y, 2) of the 
holograms at a vertical distance Zis described in Eq 1 [22]: 


16x; y; ZP Ya jY Ox; y; Zi” b jU 5x y; 2’ PU %Ox; y, ZU Ox; y; ZP p U 6x; y; ZU %x; 
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The first and second terms in Eq 1 represented the background and 
scattered light intensity. The scattered intensity is generally weaker than 
that of the background. The third and fourth terms represented the 
interference maxima and minima [22]. The CMOS sensor was placed at a 
distance Z from the light source (Z > z,) and recorded the intensity of 
the holograms [20,22]. However, the phase information was lost, which 
made it difficult to discern objects 
from a digital holographic image. 

Digital reconstruction algorithms were used to convert the holographic 
image to a micro- scope-like image. The reconstruction commonly known 
as ‘holographic reconstruction’ con- verted the holographic image in the 
object plane to the detector plane. For our application, we have used the 
Fresnel diffraction method which utilized a single Fourier transform (F) 


to back 
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Fig 4. The algorithm for the lensless SCN egg counting method is depicted as a flow chart. 
https://doi.org/10.137 1/journal.pone.0223386.g004 


propagate the wave a distance z = Z -z, to the reconstruction plane, 
UAX, Y, 2). Eq 2 describes how the reconstruction plane was obtained by 
Fourier transforming the element-wise multi- plication of image intensity 
and combining it with a Fresnel approximated transfer function and a 
phase term. This operation can be written as [23]: 
= iPx 2 ‘ Px 2 ° 

U w y;zP “= PO” Fyxy eo byr 

pp @ 5 »r®@ - 

" Iz 


where A is the wavelength and k %4 *° is the wave number. The pixel size 
An in the reconstruc- tion plane was directly proportional to the 
reconstruction plane as shown in Eq 3. 


Az 
DZ, ND 63P 
p 


where A, is the pixel pitch of the CMOS sensor and N is the minimum 
dimension of the holo- graphic image size. The reconstructed images at 
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this stage often looked blurred because of the presence of zeroth order 
frequencies, which can be removed if the phase is known. An 
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approximate phase was obtained using the Fresnel approximation. There 
are other techniques to recover the phase that involve taking multiple 
images of the same object by varying parame- ters such as wavelengths, 
angles of the LED, and z values [24,25]. Some other techniques use 
iterative methods to retrieve the phase. These alternatives are also 
effective but require addi- tional imaging components and computations. 


Results 
OptiPrep™-based density gradient centrifugation 


A qualitative analysis of the debris distribution was performed to 
characterize the effects of OptiPrep” during the centrifugation step. 
Using the protocol described in Fig 1A, different OptiPrep” solutions (0%, 
20%, 40%, 50%, 60%, and 80% by volume) were used. Images of the 
centrifuge tubes before and after the centrifugation step were taken for 
the different OptiPrep" solutions and cropped as shown in Fig 5 (i,ii). The 
three layers were separated after centrifuga- tion, diluted with water to 
12 mL, and put on separate filter papers. Fig 5 (iii) shows a plot of the 
mean pixel intensity levels as a function of distance along the center axis 
of the tube. In these plots, the orange and purple lines denote the mean 
pixel intensity levels before and after the centrifugation step, respectively. 

Compared to the control treatment (i.e. 0% OptiPrep”) (Fig 5 Control), 
the sample cleaning step using OptiPrep” resulted in a marked visual 
separation of the debris. Before centrifuga- tion, the tube had a relatively 
consistent debris distribution (Fig 5 Control i.) whereas after 
centrifugation, all of the suspended debris were concentrated into the 
pellet at the bottom of the tube (Fig 5 Control ii.). Since the control 
sample did not contain the OptiPrep” solution, there was no distinct 
interface layer. Here the mean intensity plots before and after 
centrifuga- tion are similar because no OptiPrep" solution was present 
(Fig 5 Control iii.). 

With 20% OptiPrep", there was a clear contrast between the images of 

the tube before (Fig 
5 20% i.) and after centrifugation (Fig 5 20% ii.). The mean pixel 
intensity plot showed a defini- tive change at the interface before and 
after centrifugation (Fig 5 20% iii.). The interface layer for 40% and 50% 
OptiPrep" appeared visually similar to that of the 20% concentration. 
How- ever, with increasing OptiPrep” concentrations (60% and 80%), 
there was a noticeably increas- ing amount of debris trapped at the 
interface layer both before and after centrifugation (Fig 5 60% i.-ii. and 
Fig 5 80% i.-ii.). This separation of debris at the interface layer was 
consistent with the mean color intensity plot (Fig 5 80% iii.). 


Scanner-based method: Egg recovery in the interface 
layer 


Each image of the filter paper corresponding to a different OptiPrep” 
solutions (20%, 40%, 50%, 60%, and 80%) was run through a custom 
software program, written in Python, to iden- tify the nematode eggs and 
give a count of the eggs within the image. The program initially takes a 
full image of the filter paper and then partitions it into patch images and 
sub-patches to identify and count the eggs. Fig 6A and 6B show 
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representative images of two filter papers hav- ing an uncleaned sample 
and a cleaned sample, respectively. In both cases, the program was able 
to differentiate stained eggs from the debris. 

For a given volume of the processed sample, the total SCN egg count 
was obtained by com- bining the egg counts for the top, interface, and 
bottom layers. The egg recovery ratio was cal- culated as a ratio of the 
egg count in each of the three layers to the total egg count as depicted in 
Fig 7. This was repeated for four processed samples (n = 4). The average 
egg recovery ratio and standard deviation were plotted for the top, 
interface, and bottom layers corresponding to the different OptiPrep” 
solutions. 
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Fig 5. Distribution of debris in the OptiPrep™ solution before and after centrifugation. (i., ii.) Cropped images 
of the centrifuge tube before and after centrifugation. iii.) The mean pixel intensity of the samples in the centrifuge 
tube as a function of vertical distance (in mm) was plotted for before (orange) and after (purple) centrifugation. The 
different OptiPrep” solutions used were 0%, 20%, 40%, 50%, 60%, and 80% by volume. The control (i.e. 0% OptiPrep”) 
does not have an interface layer, and so the supernatant is treated as the top and the pellet as the bottom layer. All test 


tubes pictured were 15 mL tubes. 
https://doi.org/10.1371/journal.pone.0223386.g005 


The graph in Fig 7A shows that the egg recovery ratio in the top layer 
was negligible after the OptiPrep” sample cleaning as virtually no eggs 
were found in this layer. In the interface layer, the egg recovery ratio 
improved with increasing OptiPrep” concentration and approached 95% 
for the case of the 80% concentration of OptiPrep". Meanwhile, the egg 
recovery in the bottom layer decreased with increasing OptiPrep"” 
concentrations and approached 4% for the case of the 80% OptiPrep” 
solution. For the remaining studies with OptiPrep”, the 50% OptiPrep” 
solution was used as the egg recovery at the interface layer was greater 
than 80% and the sample purity was the highest. Beyond an OptiPrep" 
concentration of 50%, sample purity started to decrease as was 
demonstrated by the increase in material remaining at the interface layer 
after cleaning (Fig 5 60% ii.). 

The graph in Fig 7A shows that most eggs were recovered from the 
interface layer of the cen- trifuge tube after sample cleaning with 
OptiPrep”. Fig 7B tabulates the average percentage of egg recovery at the 
interface layer for the associated volume percent OptiPrep” solutions (n = 
4). The spatial distribution of eggs in the images of the filter papers 
(corresponding to the interface layer) are shown as raster plots in Fig 7C. 
Four extraction runs were conducted, each with a fresh batch of extracted 
and stained SCN eggs and a pre-specified OptiPrep” solution. The raster 
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plots show the distribution of the SCN eggs (denoted as dark pixels) 
across the 324 sections of the images generated by the software. The 
raster plots show that the eggs (or dark pixels) were distributed 
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49 a 
Filter paper 


Juvenile 
from egg 


Fig 6. Software program to automatically detect SCN eggs among the debris in the scanner-based method. The software 
program starts with the full image of the filter paper and partitions it into patch images and sub-patches to identify and count the 
SCN eggs. a) An example of identifying eggs in a debris-laden, uncleaned sample is shown. b) An example of identifying eggs in an 
OptiPrep" cleaned sample is shown where it is even possible to differentiate between stained eggs and juvenile nematodes which 


have recently emerged from the eggs. Scale bar = 15 mm. 
https://doi.org/10.137 1/journal.pone.0223386.g006 
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Fig 7. Sample cleaning using OptiPrep™ and SCN egg recovery from the three layers. a) The graph shows the ratio of SCN 
eggs recovered from the three layers (top, interface, and bottom layers) at different concentrations (volume percentages) of 
OptiPrep". At 50% and higher OptiPrep"™ concentrations, the egg recovery ratio was greater than 80% in the interface layer. b) The 
egg recovery percentage in the interface layer is listed for different percentages of OptiPrep” by volume for four separate runs of 
processed samples. c) Raster plots help to visualize the spatial distribution of the eggs on the filter papers (corresponding to the 
interface layer) across four extraction runs as the volume percentage of OptiPrep" was varied. The egg count for each patch image 


is denoted as a dark pixel in the raster plot. 


https://doi.org/10.137 1/journal.pone.0223386.g007 


across most of the filter paper, which helped facilitate the identification of 
individual eggs. Com- paring the raster plots for 20% OptiPrep” with 
those of higher OptiPrep” concentrations, the ras- ter plots appear darker 
as the OptiPrep” concentration increases, suggesting that the egg 
recovery ratio increased in the interface layer as OptiPrep” 
concentrations increased. 


Performance of sample cleaning methods: OptiPrep™ 
versus sucrose 


The OptiPrep” sample cleaning method was compared with the sucrose 
centrifugation tech- nique, which is the conventional method of separating 
nematodes and eggs from debris [26]. The samples were treated as 
described in Fig 8(A). A stock sucrose solution was made by dis- solving 
454 gm of sucrose in 1 L of tap water. Approximately 500 cm? of SCN 
infested soil was manually processed (using techniques in the ‘Sample 
Preparation’ section) to extract the SCN cysts, then eggs and debris. The 
samples (i.e. liquid containing eggs and debris) were diluted to 500 mL and 
mixed thoroughly before obtaining the manual egg counts for this bulk 
sample using a nematode counting slide and dissecting microscope with a 
magnification of 100~x. 
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The bulk egg sample was stained with acid fuchsin and split into ten 
different 50 mL test 

tubes. The contents in five of these test tubes were cleaned using the 
OptiPrep” cleaning method (using a 50% by volume OptiPrep" solution). 
The egg counts for these five samples were obtained using the lensless 
method and the scanner-based egg counting method. The egg counts were 
reported as the number of eggs per milliliter. 
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Fig 8. Comparison of the sample cleaning and SCN egg counting methods. a) i-iii) Steps performed to extract cysts, then eggs 
from cysts, stain eggs, perform sample cleaning (OptiPrep" and sucrose cleaning methods), and egg counting (manual, scanner- 
based, and lensless imaging). b) The table lists the average egg count (n = 5) and standard deviation for the egg count from scanner- 


based and lensless methods compared to the standard nematode slide count on the uncleaned bulk sample. 
https://doi.org/10.137 1/journal.pone.0223386.g008 


The remaining five test tubes of samples were cleaned by the sucrose 
centrifugation tech- nique. Each tube was centrifuged for 5 minutes at 
420 G and the excess water was decanted off. The pellet was thoroughly 
mixed in 50 mL of the stock sucrose solution and centrifuged for one 
additional minute to separate the eggs from the soil pellet. The sucrose 
solution superna- tant, containing the eggs, was poured over a 37-yum- 
pore sieve and rinsed in water to remove the excess sucrose. Three 1 mL 
sub-samples of the processed egg samples were placed on a nematode 
counting slide, and the SCN eggs were manually counted under a 
microscope and averaged to estimate the egg count per milliliter. 

Soil samples were collected from two fields in Iowa-one each in 
Muscatine and Story Counties. Soil was collected from two different 
places in Iowa to get a range of soil textures with which to test the egg 
cleaning methods. Both methods of sample cleaning (i.e. OptiPrep” and 
sucrose based) were conducted on ten different samples (five from 
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Muscatine and five from Story). The egg counting was done by three 
methods-nematode counting slide, scanner- 
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based, and lensless based. The mean and standard deviation of the egg 
count obtained from the different sample cleaning methods (i.e. bulk, 
sucrose, OptiPrep”) and counting methods (i.e. nematode counting slide 
(Chalex, LLC, Park City, UT), scanner, lensless setups) are tabu- lated in 
Fig 8(B). 


Lensless imaging method: Reconstruction of holographic 
videos 


After sample cleaning using 50% OptiPrep” by volume, a 1 mL solution 
of the interface layer was loaded into a syringe. The dispensing needle of 
the syringe was connected to the input port of the microfluidic flow chip. 
The syringe pump was turned on at a flow rate of 1 mL per hour. As the 
liquid and particles began to pass through the flow channel, the RPi was 
turned on and the holographic video was recorded (Si Video). At the end 
of the video recording, the video was transferred to a remote workstation 
and processed by a MATLAB script to generate the reconstructed video. 
Reconstruction of a single frame of the video took approximately 4 
seconds on the CPU (Intel Xeon E5, 32GB RAM). 

Fig 9A illustrates the enhancement in object clarity after image 
reconstruction. Fig 9A(i.) shows a section of the image captured from the 
lensless imaging setup (i.e. raw image). After reconstruction, the clarity 
and focus of the raw image was considerably enhanced, and the 
reconstructed image is shown in Fig 9A(ii.). For comparison, the bright- 
field image of the same area was captured by a stereo microscope at 50x 
magnification and is shown in Fig 9A(iii.). 

The reconstructed video had sufficient resolution and clarity to identify 
the nematode eggs in every frame. Fig 9B shows time-lapsed image 
frames of a representative video that was recorded from the lensless 
imaging setup and reconstructed thereafter. The image frames refer to a 
small area of the flow channel. The eggs detected in each frame are 
marked in cyan color with their distinct egg number. The video was 
recorded at 1 frame per second and the flow rate was set to 1 mL/hour 
(S2 Video). 


Performance of automated egg counting methods 


We compared the accuracy of the egg counts obtained from the software 
program with manual egg counts as shown in Fig 10. In total, there were 
fifteen images of filter papers corresponding to the five different 
OptiPrep" concentrations and their three individual layers (i.e. top, inter- 
face, and bottom). Initially, the images of filter papers were divided into 
324 patch images (1024 x 1024 pixels) and saved in separate folders on 
the computer. To obtain the manual egg counts, a user viewed each 
patch image and counted the number of eggs (Fig 8). This process was 
repeated for all fifteen of the filter paper images. Next, to obtain egg 
counts from the soft- ware program, each patch image was further split 
into 256 x 256 pixel sub-patches and ana- lyzed using the trained deep 
learning model. Since the deep learning model was trained for 128 x 128 
pixel images, the sub-patches were resized before being passed to the 
learning model. Fig 10A shows the total egg counts calculated by the 
software program and by the man- ual method for all the filter paper 
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images. The R’ value between these two sets of data was high indicating 
that there was good agreement between the two sets of counts. We also 
investigated the correlation between the software and manual egg counts 
at a specific OptiPrep” concentra- tion (50%). There were a total of 972 
patch images across the three layers, and the correlation between them 
was high as well, as shown in Fig 10B. This high correlation 
demonstrates that our software program can provide the egg count in 
processed samples with accuracy that is comparable to visual detection. 
To test the correlation of the software for egg counting used with the 
lensless imaging setup, multiple samples (n = 12) were used. The 
processed samples, after cleaning with 50% OptiPrep", were passed 
through the lensless imaging setup. The raw holographic videos were 
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Fig 9. Reconstruction of the holographic videos and a time-lapsed visualization of the lensless imaging 
method. a) i) The raw image of a sample area was recorded by the lensless imaging setup (z = 8 mm). ii) The 
reconstructed image was produced using the Fresnel diffraction method. iii) The sample area was imaged with a 
stereo microscope with bright-field illumination at 50x magnification. Scale bar = 100 um. b) The images were taken 
from a representative holographic video recorded by the lensless imaging setup and reconstructed thereafter. A small 
area of the flow channel is shown here. The eggs detected by the software were labeled with a distinct egg number. 
Scale bar = 200 um. 


https://doi.org/10.1371/journal.pone.0223386.g009 


recorded as described earlier. After reconstruction, the videos were run 
through the software program to identify the eggs in every frame and 
provide the total egg count for the entire video. The processed sample 
was retrieved at the outlet of the flow channel and passed through the 
scanner-based egg counting method. The egg counts from the lensless 
method and the scanner-based method were plotted for all the samples 
as shown in Fig 10c. There was high correlation between the egg counts 
from both the methods with egg counts determined by direct microscopic 
observation. 
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Fig 10. Performance of the software programs for automated SCN egg counting. a) The plot shows the SCN egg counts 
obtained for the three layers (top, interface, and bottom) for all different OptiPrep" concentrations. The data points circled in 
orange belong to samples with fewer eggs. Similarly, the points circled in blue belong to samples with a large number of eggs. b) 
The plot shows the correlation between manual and software egg counts for all image patches of the three layers at 50% OptiPrep" 
by volume. The data were collected from three filter papers, each subdivided into 324 patch images. c) The egg counts obtained 
from the two methods are compared: software program for lensless imaging method and the scanner-based egg counting method. 
Each processed sample was passed through the lensless imaging setup and further analyzed by the scanner-based method. There is 


PLOS ONE | https://doi.org/10.1371/journal.pone.0223386 October 31/ 


‘O-PLOS ION E Removing soil debris and high-throughput nematode egg 


high correlation between the egg counts from the two methods. The plotted data corresponds to 12 different processed samples. 
https://doi.org/10.1371/journal.pone.0223386.g010 
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Discussion 
OptiPrep™-based sample cleaning method 


The current state-of-the-art sample cleaning method involves the use of 
sucrose centrifugation to separate plant-parasitic nematode juveniles and 
eggs from debris. While sucrose is inexpen- sive and readily available, its 
high osmotic potential may have adverse effects on nematode juveniles 
and eggs upon prolonged exposure. Also, if the sucrose solution is not 
adequately washed away from the nematodes and eggs, there is the risk 
of bacterial and fungal growth in the samples. To address these 
challenges of sucrose centrifugation, we developed an alternative density 
gradient method using OptiPrep™ solution to separate plant-parasitic 
nematode eggs from debris. OptiPrep™ has certain advantages over 
sucrose. OptiPrep™ is iso-osmotic and nontoxic to biological materials 
such as cells and tissues. OptiPrep™ has a lower viscosity compared to 
sucrose, which eases the mixing and washing steps with water. The 
lower viscos- ity and osmolality of OptiPrep™ is advantageous to clean 
small volumes of sample homoge- nates (5-10 mL in a 15 mL test tube) 
and to condense the eggs into the interface layer (2 mL) within the 
centrifuge tube. Greater nematode egg recovery is achieved from 40% to 
80% Opti- 

Prep™ solutions (i.e. > 80% recovery) compared to those from the sucrose 
centrifugation (i.e. 

less than 40% recovery) (Fig 7). As a disadvantage, the price of the 
OptiPrep™ density gradi- ent medium is higher (approximately $250 USD 
for 250 mL) than sucrose. The high price of OptiPrep™ may be a 
challenge for soil processing facilities which perform large-scale or rou- 
tine separation of nematode eggs from debris. 


Automated nematode egg counting methods 


The traditional method of nematode egg counting is based on 
microscopic observation of sam- ples. The chamber of nematode counting 
slides is loaded with 1 mL of the sample, and the number of eggs on the 
grid is counted manually. Counting nematode eggs in this manner is a 
cheap and relatively simple procedure but requires considerable time, 
expertise, and patience of a trained expert [26]. Also, cross 
contamination can occur because of insufficient or improper cleaning of 
the chamber in the nematode counting slide between samples. A deep 
learning architecture was recently developed to identify SCN eggs from 
microscopic images of nematode-counting slides [27]. Their input images 
(480 x 640 pixels) were subdivided into 
16 x 16 image patches and passed through a selectivity function to 
identify patches with an 
SCN egg. However, with this deep learning architecture, it is challenging 
to handle situations with larger image sizes or multiple eggs (i.e. greater 
than two eggs) within an image patch such as those obtained from bulk 
samples on our high-resolution scanner. 

The new, scanner-based egg counting method we developed has 
advantages over previous counting methods. Here the filter papers are 
automatically scanned to obtain the egg counts using a trained deep 
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learning model. Our deep learning model handles large input images 
(17759 x 17759 pixels) which are subdivided into 4900 patches (256 x 
256 pixels). The U-Net convolutional autoencoder model can recognize 
multiple eggs within a patch using object localization. The time to obtain 
the computerized egg counts from a single filter paper image is around 5 
minutes, whereas this time is much longer for manual counting and varies 
across dif- ferent samples and persons. While the size of the sample 
assessed using a nematode counting slide is limited by the volume of the 
observation chamber (1 mL) of the slide, each filter paper can process 
around 12 mL of the sample suspension. Also, cross contamination 
between sam- ples is eliminated in the scanner-based method as the 
substrate is a disposable filter paper. The filter paper substrate is 
relatively inexpensive (approximately $0.48 USD) compared to the 
nematode counting slide (approximately $50 USD). As a disadvantage, the 
scanner-based 
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method requires a pre-processing step where the sample needs to be 
stained for proper identi- fication. The staining process is performed at 
high temperature which eventually kill the eggs. Also, it is difficult to 
recover the eggs after imaging as the processed sample is distributed on a 
filter paper. Hence the scanner-based method is not suited for 
experiments where samples containing live eggs or nematodes are to be 
counted for further use, such as to infest soil for greenhouse experiments. 
The new, lensless egg counting method has advantages over the 
manual slide counting and scanner-based methods. Here egg staining is 
not required as the lensless imaging setup can enumerate both stained 
and unstained samples. Because of continuous fluid flow through the 
ports, the counted eggs are easily recovered for later use from the 
microfluidic flow chip. The material cost is relatively low (less than $100 
USD) for the lensless imaging setup and adhesive tape microfluidics 
compared to the costs of standard microscopes, scanners, and nematode 
counting slides used in previous methods. The lensless imaging method 
is functionally more automated than the scanner-based method as a user 
only needs to load a syringe containing the egg sample into a syringe 
pump, set the right volume and speed parameters, and start the video 
recording process. As a disadvantage, the lensless egg counting method 
has a relatively long processing time. The flow rate of the sample is 1 mL 
per hour; higher flow rates produce noise in the recorded images. The 
reconstruction of the holographic video takes a few hours and the egg 
counting software needs another 10 minutes to obtain the egg count. 


Conclusion 


In conclusion, we developed methods to improve the efficiency of SCN 
egg counting by sam- ple cleaning and imaging on scanner and lensless 
setups. We also created software programs to count the eggs. The use of 
OptiPrep" during density gradient centrifugation helped capture most of 
the nematode eggs in the interface layer. Egg recovery was greater than 
80% in the interface layer for the case of 50% OptiPrep” by volume. 
Then the processed samples (i.e. eggs with debris) were either placed on 
filter papers to record static images with a flatbed scanner or passed 
through a microfluidic flow chip to record real-time, holographic videos 
by a lensless imaging setup. In the scanner-based method, a 
convolutional autoencoder network recognized the nematode eggs from 
static scanner images with reasonable accuracy as confirmed by visual 
observations. For the lensless method, the software program successfully 
reconstructed the holographic videos and identified the nematode eggs in 
the processed sample, which was sub- sequently confirmed by the 
scanner-based method. The cost of the materials was low by using off- 
the-shelf components (such as CMOS sensor, Raspberry Pi, LED), 
microfluidic flow chips made from inexpensive double-sided adhesive 
tapes, and standard filter papers. Furthermore, the use of image 
processing and deep learning tools circumvents the need for hiring and 
train- ing personnel to count nematode eggs using nematode counting 
slides. With our new meth- ods, manual intervention is only needed 
during sample preparation and loading; the software programs handle 
the remaining operations automatically (i.e. image or video capture, 
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detec- tion and counting, data storage) over a broad range of nematode 
egg numbers (10-300 eggs/ mL). 
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